Weak ferromagnetism with the Kondo screening effect in the Kondo lattice systems 
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We carefully consider the interplay between ferromagnetism and the Kondo screening effect in the 
conventional Kondo lattice systems at finite temperatures. Within an effective mean-field theory 
for small conduction electron densities, a complete phase diagram has been determined. In the 
ferromagnetic ordered phase, there is a characteristic temperature scale to indicate the presence of 
the Kondo screening effect. We further find two distinct ferromagnetic long-range ordered phases 
coexisting with the Kondo screening effect: a spin polarized phase and a non-polarized phase. A 
continuous phase transition exists to separate the spin non-polarized ferromagnetic ordered phase 
from the paramagnetic heavy Fermi liquid phase. These results may be used to explain the weak 
ferromagnetism observed recently in the Kondo lattice materials. 
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The most important issue in the study of heavy fermion 
materials is the interplay between the Kondo screen- 
ing and the magnetic interactions among local mag- 
netic moments mediated by the conduction electrons^— 
The former effect favors the formation of Kondo sin- 
glet state in the strong Kondo coupling limit, while 
the latter interactions tend to stabilize a magnetically 
long-range ordered state in the weak Kondo coupling 
limit. In-bctwccn these two distinct phases, there ex- 
ists a magnetic phase transition. Although such a phase 
transition was suggested by Doniach many years ago^i 
the complete finite temperature phase diagram for the 
Kondo lattice systems has not been derived from a mi- 
croscopic theory^ At the half-filling of the conduction 
electrons, the antiferromagnetic long-range order dom- 
inates over the local magnetic moments, which can be 
partially screened by the conduction electrons in the 
intermediate Kondo coupling regime* 7 - - — Very recently, 
close to the magnetic phase transition, weak ferromag- 
netism below the Kondo temperature has been discov- 
ered in the Kondo lattice materials UCub x Pd, T (ReffTlh. 
URhi^RuzGe (RefGJ),YbNi 4 P 2 (Ref[m YbCu 2 Si 2 
(ReflU), and Yb(Rh .73Coo.27)2Si 2 (Refill). So an in- 
teresting question arises as whether the ferromagnetic 
long-range order can coexist with the Kondo screening 
effect. 

To account for the ferromagnetism within the Kondo 
lattice model, one can assume that conduction electrons 
per local moment n c is far away from half-filling, where 
the ferromagnetic correlations dominate in the small 

Kondo coupling regime^ia In our previous paper^ we where £ k is the dispersion of the conduction electrons, 



resulting ground state is either a spin-polarized or a non- 
polarized ferromagnetic phase. The existence of the spin- 
polarized phase has been confirmed by the recent dy- 
namic mean-field calculations in infinite dimensions and 
density matrix renormalization group in one dimension, 
where such a state is referred to as the spin-selective 
Kondo insulator^ 

In this paper, we will further explore the global phase 
diagram of the Kondo lattice model at finite tempera- 
tures for small conduction electron densities. Below the 
Curie temperature, we find a characteristic temperature 
scale to signal the Kondo screening effect for the first 
time. Moreover, there exist a spin polarized phase and 
a non-polarized ferromagnetic long-range ordered phase 
coexisting with the Kondo screening effect. The former 
phase spans a large area in the phase diagram, while the 
latter phase just sits in a very narrow region close to 
the phase boundary with the paramagnetic heavy Fermi 
liquid phase. Moreover, a second order phase transition 
occurs from the spin non-polarized ferromagnetic ordered 
state to the paramagnetic heavy Fermi liquid state, and 
the transition line becomes very steep close to the quan- 
tum critical point. Our results may be used to explain 
the weak ferromagnetism and quantum critical behavior 
observed in YbNi4P 2 

The Hamiltonian of the Kondo lattice systems is de- 
fined by 



H = ^2 eucl a c ka + J K ^2 °i ' S j 



(1) 
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have carefully studied the possible ground state phases 
within an effective mean field theory. In particular, for 
0.16 < n c < 0.82 and close to the magnetic phase tran- 
sition, the local moments can be only partially screened 
by the conduction electrons, and the remaining uncom- 
pensated parts develop the ferromagnetic long-range or- 
der. Depending on the Kondo coupling strength, the 



Cj = | J2a0 c la T aP c iP is the spin density operator of 
the conduction electrons, r is the Pauli matrix, and 
the Kondo coupling strength Jjc > 0. When the lo- 



calized spins are denoted by S, = ^J2 a p fL T apL 
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the pseudo-fermion representation, the projection onto 
the physical subspace has to be implemented by a local 
constraint ^ Haf™ = 1- ^ is straightforward to dc- 
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FIG. 1: (Color online) The schematic phase diagram expected 
from the interplay between ferromagnetic correlations and the 
Kondo screening effect. T<£ denotes the Curie temperature in 
the absence of the Kondo effect, and T% represents the Kondo 
temperature without the ferromagnetic correlations. 



(cL/j^ + fi±Cii) is introduced to decouple the transverse 
exchange term. We also introduce a Lagrangian multi- 
plier A to enforce the local constraint, which becomes 
the chemical potential in the mean field approximation. 
Then the mean field Hamiltonian in the momentum space 
can be written in a compact form 



(2) 

where e kff = e k + ^pa, X a = X + ±^a, e Q = 

JK"m c mj — A, a = ±1 denote the up and down 

spin orientations, and M is the total number of lattice 
sites. The quasiparticlc excitation spectra are thus ob- 
tained by 



compose the Kondo spin exchange into longitudinal and 
transversal parts 

a* ■ S, = a!Sf - \i( C y it + fl iCli f + (cy n + f} fCif n 

where the longitudinal part describes the polarization of 
the conduction electrons, giving rise to the usual RKKY 
interaction among the local moments; while the trans- 
verse part represents the spin-flip scattering of the con- 
duction electrons by the local moments, yielding the local 
Kondo screening effecti^ The competition between these 
two interactions determines the possible ground states of 
the Kondo lattice systems. 

Similar to the interplay between the antiferromagnctic 
correlations and the Kondo screening effect argued by 
Doniach,— a schematic finite temperature phase diagram 
can be expected for the interplay between the ferromag- 
netic correlations and the Kondo screening effect. In 
Fig.l, the Curie temperature is plotted as a function 
of the Kondo coupling. For the small Kondo couplings, 
the ferromagnetic ordering (Curie) temperature is larger 
than the single-impurity Kondo temperature. When the 
Kondo coupling is strong enough, the Curie temperature 
is suppressed completely. However, there is an important 
issue as whether there should be a characteristic temper- 
ature scale inside the ferromagnetic ordered phase to sig- 
nal the presence of the Kondo screening effect. If so, there 
may exist two distinct ferromagnetic ordered phases: a 
pure ferromagnetic phase with a small Fermi surface con- 
sisting of conduction electrons only, and a ferromagnetic 
phase with an enlarged Fermi surface including both con- 
duction electrons and local magnetic moments, coexisting 
with the Kondo screening. 

Let us first review the effective mean field theory 
for the ground state used in our previous study— We 
introduce two ferromagnetic order parameters: mj = 
(Sf) and m c — (erf) to decouple the longitudinal ex- 
change term, while a hybridization order parameter V = 



eua + A CT ± J (e kff - A(j ) + ( J K V) 2 



(3) 



where there appear four quasiparticlc bands with spin 
splitting. 

Using the method of equation of motion, the single 
particle Green functions can be derived, while the cor- 
responding density of states can be calculated and ex- 
pressed as 



p a c (uj) = p° c [9(uj - uj 1<7 )9(lo 2 v - u) + 6{uj - uj 3a )9(uj 4a - u)], 



where 0{oj) is a step function and a constant density of 
states of conduction electrons has been assumed p° c = , 
with D as a half- width of the conduction electron band. 
The four quasiparticlc band edges can be expressed as 



^Icr = 

^2cr = 

^3<7 = 

W4„ = 

where 
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Then using the spectral representation of the Green 
functions, we derive the mean- field equations at finite 
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FIG. 2: Schematic DOS in the presence of Kondo screening 
effect, (a) for the spin non-polarized ferromagnetic phase, (b) 
for the spin polarized ferromagnetic phase. 



temperatures as follows 



[pt(u) + p c (w)] = n c 
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ordered phases with a finite value of the hybridization 
parameter V may appear, depending on the Kondo cou- 
pling strength. We call them polar-II and non-polar-II 
ferromagnetic ordered states, respectively. For a strong 
Kondo coupling, the pure Kondo paramagnetic phase 
is the ground state. There is a continuous transition 
from the non-polar-II ferromagnetic ordered phase to the 
Kondo paramagnetic phase. 

Now we calculate the finite temperature phase dia- 
gram. First of all, if the temperature is high enough, all 
order parameters must disappear, so the conduction elec- 
trons and local moments are decoupled. As the temper- 
ature is decreased down to the Curie temperature of the 
pure ferromagnetic phase TjS, both m c and m/ approach 
zero, but the ratio m c /mf is finite. The self-consistent 
equations give rise to 



x ~ Jk 



-g^r m c ( 6 ) 



and the Curie temperature Tq can be estimated as 

J 2 

J-n — 



32D' 



(7) 



which is independent of the density of conduction elec- 
trons. 

On the other hand, if Jk is large enough, the sys- 
tem must be in pure Kondo paramagnetic phase. As the 
Kondo coupling decreases, the Kondo singlets are desta- 
bilized. When T — > T K , the hybridization vanishes, and 
the self-consistent equations can be reduced to 



where f(uj) — 1/ [l + e^ u ~ t ^' T ] is the Fermi distribu- 
tion function. To make the magnetic interaction be- 
tween the nearest neighboring local moments ferromag- 
netic, we should confine the density of conduction elec- 
trons to n c < 0.82 from the previous mean field study.— 

The position of the chemical potential p with respect to 
the band edges is very important. At zero temperature, 
there are two different situations. The corresponding 
schematic local density of states are displayed in Fig. 2. 
For wi_ < p < 0J2+, both the lower spin- up and spin- 
down quasiparticle bands are partially occupied, corre- 
sponding to the spin non-polarized state. However, for 
UJ2+ < p < 0J2—, the lower spin-up quasiparticle band 
is completely occupied, while the lower spin-down quasi- 
particle band is only partially occupied, corresponding 
to the spin polarized stated An energy gap A-f exists in 
the spin-up quasiparticle band, and there is a plateau in 
the total magnetization: m c + m/ = (1 — n c )/2. 

The ground-state phase diagram has been obtained in 
our previous study^ When n c < 0.16, the polarized pure 
ferromagnetic (polar-I) phase is a ground state in the 
large Kondo coupling region. For 0.16 < n c < 0.82, 
the ground state is given by the non-polarized pure fer- 
romagnetic (non-polar-I) phase in the weak Kondo cou- 
pling limit; while in the intermediate Kondo coupling 
regime, both polarized and non-polarized ferromagnetic 
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(8) 



which can be numerically solved, leading to the Kondo 
temperature T K . 

After obtaining XjS and T K , we expect that the pure 
ferromagnetic phase exists for > T K in the small 
Kondo coupling limit, while for T K > XjS the Kondo 
screening is present, and the competition between the 
ferromagnetic correlations and Kondo screening effect 
should be taken into account more carefully. 

In the presence of the Kondo screening, the corre- 
sponding Curie temperature Tq can be still defined. As 
T — > Tq, the magnetic moments m c and m/ approach 
zero, but their ratio is finite, m c /mf ^ 0. The self- 
consistent equations Eq.(5) can be solved numerically, 
leading to the Curie temperature Tq and the mean-field 
parameters of p, A, V, and m c /mf. On the other hand, 
when the ferromagnetism is present, we can also intro- 
duce the Kondo temperature Tk incorporating the hy- 
bridization effect. When V — > and T — s- Tk, the nu- 
merical solution of the self-consistent equations gives rise 
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FIG. 3: (Color online) The finite temperature phase diagram 
at n c = 0.2. In addition to the pure ferromagnetic ordered 
phase (V — 0, m c 7^ and m/ 7^ 0) and the Kondo param- 
agnetic phase (V 7^ 0, m c = mj = 0), there are two differ- 
ent ferromagnetic ordered phases coexisiting with the Kondo 
screening (V 7^ 0, m c 7^ and mf 7^ 0): the spin polarized 
phase (A-)- 7^ 0) and spin non-polarized phase (A-f = 0). 



to the Kondo temperature Tk and the mean-field param- 
eters fi, A, m c , and m/. 

The resulting phase diagram is shown in Fig. 3 for 
n c = 0.2. As the Kondo coupling Jk increases from a 
small value, the Curie temperature Tq first increases up 
to a maximal value, and then continuously decreases to 
zero at J£? = 1.133D. For small values of Jk/D < 0.41, 
the Kondo temperature Tk vanishes. However, when 
Jk/D > 0.41, the Kondo temperature curve consists of 
two parts, meeting each other precisely at the Curie tem- 
perature (Tk = 2c). Inside the ferromagnetic ordered 
phase, Tk starts from a finite value and then decreases 
down to zero at J^- 1 = 0.5061?; while in the paramag- 
netic phase, Tk follows the behavior of the bare Kondo 
temperature T K . 

In the coexistence region, depending on whether the 
chemical potential /1 is inside the energy gap of spin-up 
quasi-particle (as shown in Fig. 2), we can calculate the 
lowest excitation energy defined by A-j- = /1 — ll>2+- In 
Fig. 4, we show Af as a function of T with fixed Kondo 
coupling parameters Jk/D = 0.6, 0.7, 0.733, 0.9 and 1.0, 
respectively. It is clearly demonstrated that the spin gap 
A-j- has a non-monotonic behavior as the temperature in- 
creases. Notice that Jk/D = 0.733 corresponds to the 
critical value between the spin polarized phase and non- 
polarized phase at zero temperature. When A-j- — > 0, the 



characteristic temperature Tq is determined, leading to 
the phase boundary separating the spin polarized and the 
non-polarized ferromagnetic ordered phases. The spin 
polarized ferromagnetic order phase spans a large area 
in the coexistence region, while the spin non-polarized 
phase just sits in a narrow strip close to the phase bound- 
ary with the paramagnetic heavy Fermi liquid phase. The 
existence of the non-polarized ferromagnetic order phase 
should be expected. However, it has been missed by the 



FIG. 4: (Color online) The energy gap of spin-up quasiparti- 
cles A-f as a function of temperature T in the coexisting phase 
at n c = 0.2. 
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FIG. 5: (Color online) The ferromagnetic magnetizations and 
hybridization parameters as a function of the Kondo coupling 
Jk with a fixed temperature at n c = 0.2. (a) T = 0.0025D, 
(b) T = 0.0075D. 



recent numerical calculations.— 

Moreover, the magnetizations of the local moments 
and the conduction electrons mf and m c are calculated 
as a function of the Kondo coupling strength Jk for 
T = 0.0025Z? and 0.0075D, as shown in Fig.5., respec- 
tively. It is clear that m c has an opposite sign of m/, 
due to the antiferromagnetic coupling between the lo- 
cal moments and conduction electrons. In order to dis- 
play the Kondo screening effect, we have also plotted the 
hybridization parameter V as a function of Jk for the 
same fixed temperatures. For low temperatures shown 
in Fig. 5a, the Kondo screening effect emerges inside the 
ferromagnetic ordered phase, and a small jump is in- 
duced in both magnetizations mf and m c . When the 
magnetization vanishes, the hybridization V has a cusp. 
However, for high temperatures shown in Fig. 5b, the fer- 
romagnetic ordering appears inside the Kondo screened 
region. The cusps in the hybridization curve are induced 
when the ferromagnetic order parameters start to emer- 
gence or vanish. 
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FIG. 6: (Color online) The magnetizations and hybridization 
parameter as a function of temperature for a given Kondo 
coupling strength at n c = 0.2. (a), (b), (c), (d), (e), and 
(f) correspond to Jk/D = 0.25, 0.45, 0.60, 1.0, 1.1, and 1.2, 
respectively. 



The magnetizations m c and m / and hybridization pa- 
rameter V have also been calculated as a function of 
temperature T for the fixed Kondo coupling strength Jk, 
which is shown in Fig. 6. For a small value of Jk / D = 0.2, 
as the temperature is increased, the magnetic moments 
m c and m/ in the absence of the Kondo screening de- 
crease down to zero at the Curie temperature Tc (see 
Fig. 6a). In contrast, for a large value of Jk/D = 1.2, 
the system is in a paramagnetic heavy Fermi liquid phase 
without ferromagnetic order (see Fig.6f). 

For Jk/D = 0.45, the Kondo screening effect starts to 
appear in the presence of ferromagnetic ordering. When 
the ferromagnetic order disappears at Tc, the hybridiza- 
tion reaches a maximal value and then decreases down to 
zero at Tk (see Fig. 6b). For larger values of the Kondo 



coupling Jk/D = 0.6, 1.0, and 1.1, the Kondo screen- 
ing effect dominates in the temperature range, and the 
ferromagnetic ordering phase occurs only in a small re- 
gion, displayed in Fig. 6c, Fig.6d and Fig.6e, respectively. 
These three figures demonstrate the interplay between 
the Kondo screening effect and the ferromagnetic corre- 
lations in the presence of thermal fluctuations. 

Finally, it is important to mention a new Kondo lat- 
tice system YbNi4P2, recently discovered by distinct 
anomalies in susceptibility, specific heat and resistivity 
measurements^ Growing out of a strongly correlated 
Kramers doublet ground state with Kondo temperature 
Tk ~ 8K, the ferromagnetic ordering temperature is 
severely reduced to T c = 0.17K with a small magnetic 
moment mj ~ 0.05/xb. Here we would like to explain 
the small ferromagnetically order moment and the sub- 
stantially reduced Curie temperature as originating from 
the presence of the Kondo screening effect, see Fig. 6c, 
Fig.6d, and Fig. 6c. The experimental results can cer- 
tainly be understood in terms of our effective mean field 
theory. The quantum critical behavior observed exper- 
imentally requires a quantum critical point separating 
the ferromagnetic ordered phase from the Kondo para- 
magnetic phase at zero temperature, which is also con- 
sistent with our finite temperature phase diagram. The 
further detailed calculations concerning with thermody- 
namic properties of the heavy fermion ferromagnetism 
arc 1 left for our future research. 

In summary, within an effective mean-field theory for 
small conduction electron densities, we have derived the 
finite temperature phase diagram. Inside the ferromag- 
netic ordered phase, a characteristic temperature scale 
has been found to signal the Kondo screening effect for 
the first time. In additional to the pure ferromagnetic 
phase, there are two distinct ferromagnetic long-range 
ordered phases coexisting with the Kondo screening ef- 
fect: a spin polarized phase and a non-polarized phase. 
A second-order phase transition and a quantum critical 
point have been found to separate the spin non-polarized 
ferromagnetic ordered phase and the paramagnetic heavy 
Fermi liquid phase. To some extent, our mean field the- 
ory has captured the main physics of the Kondo lattice 
systems, which provides an alternative interpretation of 
weak ferromagnetism observed experimentally 
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